Coulomb stability of the 47r-periodic Josephson effect of Majorana fermions 
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The Josephson energy of two superconducting islands containing Majorana fermions is a 47r- 
periodic function of the superconducting phase difference. If the islands have a small capacitance, 
their ground state energy is governed by the competition of Josephson and charging energies. We 
calculate this ground state energy in a ring geometry, as a function of the flux $ enclosed by the ring, 
and show that the dependence on the Aharonov-Bohm phase 2e$/ft remains 47r-periodic regardless 
of the ratio of charging and Josephson energies — provided that the entire ring is in a topologically 
nontrivial state. If part of the ring is topologically trivial, then the charging energy induces quantum 
phase slips that restore the usual 27r-periodicity. 

PACS numbers: 73.23.Hk, 74.50.+r, 74.78.Na, 74.81. Fa 



The energy Hj of a tunnel junction between two su- 
perconductors (a Josephson junction) depends on the dif- 
ference (j) of the phase of the order parameter on the two 
sides of the junction. The derivative Ij — {2e/h)dHj/d(f> 
gives the supercurrent flowing through the junction in 
the absence of an applied voltage. In a ring geometry, 
the supercurrent depends periodically on the flux $ en- 
closed by the ring, with periodicity h/2e. This familiar 
DC Josephson effect^'^ acquires a new twist if the junction 
contains Majorana fermions<^"— 

Majorana fermions are charge-neutral quasiparticles 
bound to midgap states, at zero excitation energy, 
which appear in a socalled topologically nontrivial 
superconductor While in the conventional Joseph- 
son effect only Cooper pairs can tunnel (with probability 
r <IC 1), Majorana fermions enable the tunneling of sin- 
gle electrons (with a larger probability ^/t). The switch 
from 2e to e as the unit of transferred charge amounts to 
a doubling of the fundamental periodicity of the Joseph- 
son energy, from Hj cx cost/) to Hj oc cos{(f>/2). In 
a ring geometry, the period of the flux dependence of 
the supercurrent Ij doubles from 2t: to An as a function 
of the Aharonov-Bohm phase* (po — 2e^/h. This 47r- 
periodic Josephson effect has been extensively studied 



theoreticallyi 
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as a way to detect the (so far, elusive) 



Majorana fermions 

Since the Majorana fermions in a typical experiment 
will be confined to superconducting islands of small ca- 
pacitance C, the Coulomb energy He = Q^/2C asso- 
ciated with a charge difference 2Q across the junction 
competes with the Josephson energy. The commuta- 
tor [(/), Q] = 2ei implies an uncertainty relation between 
charge and phase differences, so that a nonzero He in- 
troduces quantum fluctuations of cj) in the ground state.- 
What is the fate of the 47r-periodic Josephson effect? 

As we will show in this paper, the supercurrent through 
the ring remains a 47r-periodic function of ipo, regardless 
of the relative magnitude of He and Hj. This Coulomb 
stability requires that all weak links in the ring contain 
Majorana fermions. If the ring has a topologically triv- 
ial segment, then quantum phase slips restore the con- 
ventional 27r-periodicity of the Josephson effect on suffi- 
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FIG. 1: Geometry of a DC SQUID, consisting of a super- 
conducting ring (grey) interrupted by two tunnel junctions 
(black) and threaded by a magnetic flux $. A semiconduc- 
tor nanowire (yellow) contains Majorana fermions at the end 
points (red dots). The two panels distinguish the cases that 
Majorana fermions are present at both junctions (top), or 
only at a single junction (bottom). The 47r-periodic Joseph- 
son effect is stable against quantum phase slips in the first 
case, but not in the second case. 



ciently long time scales. We calculate the limiting time 
scale for the destruction of the 47r-periodic Josephson ef- 
fect by quantum phase slips and find that it can be much 
shorter than the competing time scale for the destruction 
of the 47r-periodicity by quasiparticle poisoningi^ 

We apply the general theory of Majorana- Josephson 
junction arrays of Xu and FuiS, to the DC SQUID geom- 
etry of Fig. [1] consisting of two superconducting islands 
separated by tunnel junctions. The islands have a charge 
difference 2Q — Qi — Q2, with Q„ = — 2ei9/90„ canon- 
ically conjugate to the superconducting phase 0„. The 
gauge invariant phase differences across the two junctions 
are given by (j) — (fii — (1)2 and ipo — 4>. Here we assume 
that the ring is sufficiently small that the flux generated 
by the supercurrent can be neglected, so the enclosed flux 
equals the externally applied fluxJi 

Each island contains a segment of a semiconductor 
nanowire, driven into a topologically nontrivial supercon- 
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ducting state by the proximity effect. ^'1" (Alternatively, 
the nanowire could be replaced by the conducting edge of 
a two-dimensional topological insulator j^) The Majorana 
fermions appearing at the end points of each segment 
are represented by ant i- commuting Hermitian operators 
71,72,73,74 that square to unity, 
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(1) 



The Majorana fermions are coupled by the tunnel junc- 
tion. We distinguish two cases. In the first case (top 
panel in Fig. [ij each of the two tunnel junctions couples 
a pair of Majorana fermions. In the second case (bot- 
tom panel) one pair of Majorana fermions is coupled by 
a Josephson junction, while the other pair remains iso- 
lated. 

The Hamiltonian H — He + Hj i + Hj,2 is the sum of 
charging and Josephson energies. 



He = ^(Q + gind)', 



Hj^i = £^Af,iri cos -~ Ej^i cos 0, 

Hj^2 = Em,2^2 cos Ej^2 COS(93o 

Ti = 17273, r2 = «747i- 



(2) 

(3) 

(4) 
(5) 



The induced charge qi^d — CgVg accounts for charges on 
nearby electrodes, controlled by a gate capacitance Cg 
and gate voltage Vg. The energy scales EM,n and Ej^n 
quantify the Josephson coupling strength of, respectively, 
single electrons and electron pairs. With this Hamilto- 
nian we can describe both cases considered, by putting 
Em^2 = for the junction without Majorana fermions. 

The eigenstates 5'(0i,(/)2) of H should satisfy the 
fermion parity constraints^ 



*(0i + 2^,02 + 2^™) = (-i)"''H-ir«^*(<^i,<^2), 

gn = ^(l-Pri), Pl=«7l72, P2=«7374- 



(6) 
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The operators g„ and pn have, respectively, eigenvalues 

0. 1 and ±1, depending on whether island n contains an 
even or an odd number of electrons. The constraint ^ 
enforces that the eigenvalues of Qn are even multiples of 
e for g„ = 0,p„ = 1 and odd multiples of e for ~ 

1, Pri -1- 

It is possible to solve the eigenvalue problem H'^ = 
E'i' subject to the constraint along the lines of Ref. 
[TgI . but alternatively one can work in an unrestricted 
Hilbert space. The restriction is removed by the unitary 
transformation 



^ = C/1C/2*, Un = exp(fg„0„/2). 



(8) 



The function ^'(01, is 27r-periodic in each of its ar- 
guments, so the constraint © is automatically satisfied. 
Now the eigenvalues of Qn are all even multiples of e. The 



transformed Hamiltonian H = {UiU2y HU1U2 becomes 



H 
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2CV'" ■ 2 
+ \ [e-'9i^i {Em.1^1 + SM,2r2e''^''/2)e*«^^^ + H.c] 
- Sj,! COS0 - £:j,2 cos(93o - </»), (9) 



where we have used the identity 



(10) 



Notice that the Hamiltonian has become 27r-periodic in 
the superconducting phases 0i , 02 , while remaining An- 
periodic in the flux (pQ. Notice also that H may depend 
on the 0„'s separately, not just on their difference. This 
does not violate charge conservation, because the conju- 
gate variables Qn now count only the number of Cooper 
pairs on each island — not the total number of electrons. 

The four Majorana fermions encode a qubit degree of 
freedomiiS The states of the qubit are distinguished by 
the parity of the number of electrons on each island. 
If the total number of electrons in the system is even 
{V = 1), the qubit states are |11) and |00), while for an 
odd total number of electrons {V = —1) the states are 
1 10) and |01). In this qubit basis, the products of Ma- 
jorana operators appearing in the Hamiltonian ^ are 
represented by Pauli matrices. 
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Tl = — (Tx, r2 — VUx- 

(11) 

It is straightforward to calculate the eigenvalues of 
H, by evaluating its matrix elements in the basis of 
eigenstates of Q. The spectrum E^^ {ipo , qind) as a func- 
tion of the enclosed flux and the induced charge has 
two branches distinguished by the total fermion parity 
V = ±1, with 



En{VQ,qhid) = K„ (<^o + 27r, gind + e/2). 



(12) 



We first consider the case that both junctions contain 
Majorana fermions (top panel in Fig. [T]). 

A fully analytical calculation is possible in the limit 
that the charging energy dominates over the Josephson 
energy {Ec = e^/2C > Em.u, Ej,n)- Only the two eigen- 
states of Q with lowest charging energy EzL^S are needed 
in this limit and 2e tunnel processes may be neglected rel- 
ative to e tunnel processes (so we may set Ej^n = 0). We 
thus obtain the simple expression 



EV^E±h 



-E 



MA' 



-EIi^2 + ^VEm,iEm,2 cos^ 



1/2 



(13) 

The resulting 47r-periodic spectrum is shown in Fig. [21 

The crossing of the two branches E^ and EZ at ipo — tt 
is protected, regardless of the value of Ec, because 
the charging energy cannot couple states of different V. 
Quasiparticle poisoning (the injection of unpaired elec- 
trons) switches the fermion parity on a time scale Tp, 
which means that the 47r-periodicity of the energy of the 
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FIG. 2: Spectrum of the DC SQUID in the top panel of Fig. [T] 
containing Majorana fermions at both Josephson junctions. 
The curves are the result in the limit that the charging 

energy dominates over the Josephson energy. The parameters 
chosen are -Em.i = Em.2 = S. The level crossing is between 
states of different fermion parity V, and therefore there can 
be no tunnel splitting due to the Coulomb interaction (which 
conserves V). 

ring can be observed if the enclosed flux is increased by 
a flux quantum in a time T$ <C Tp. 

We now turn to the case that one of the two Joseph- 
son junctions does not contain Majorana fermions (lower 
panel in Fig. [1]). By putting Em,2 = the Hamiltonian 
becomes 27r-periodic in ipo. In Fig. |3]we show the spec- 
trum for a relatively large Josephson energy of the trivial 
junction. The phase is then a nearly classical variable, 
which in the ground state is close to (po (mod 2tt). The 
charging energy opens a gap in the spectrum near ipQ — t: 
(mod 27r), by inducing tunnel processes from (f> = ipQ to 
4> = ipq±2tt (quantum phase slips) . A tunnel splitting by 
the P-conserving charging energy is now allowed because 
the level crossing is between states of the same V. 

A semiclassical calculation of the tunnel splitting due 
to quantum phase slips at the trivial Josephson junction, 
along the lines of Ref. HO, gives for Ej = Ej^2 ^ Ec ^ 
Em,i = Em the spectrum 

El = -Ej + ^2EcEj ±^EIj cos2(^o/2) + A^, (14) 
A = 16{EcEy2n^Y^\xp{-^8Ej/Ec) 



X y cos2(^q;„,/e) + sin^nqUe), (15) 

where we have abbreviated q[^^ — (7ind+(e/4)(l— T'). The 
second term on the right-hand-side of Eq. (|14p describes 
the effect of zero-point fluctuations of around the val- 
ues ipo and ipo ± 2tt. Tunnel processes (j) — tpo '-^ fo + 2Tr 
and (f> = ifo t-^ ipo — 2t: produce the third term. The 
sine and cosine factors in Eq. (jlSp accounts for inter- 
ference between these two quantum phase slip processes 
(Aharonov-Casher effect) i^"— The numerical calculation 
in Fig. 2] agrees quite well with the semiclassical approx- 
imation (IT5|) . 
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FIG. 3: Spectrum of the DC SQUID in the bottom panel of 
Fig. [T] containing Majorana fermions at only one of the two 
Josephson junctions. The curves are a numerical calculation 
for the full Hamiltonian, in the regime that the Josephson 
energy of the trivial junction is the largest energy scale. The 
parameters chosen are Ej^2 = 4:Ec ~ 10-EAf,ii Em,2 = = 
Ej^i, and gind ~ 0. In contrast to Fig. [2l a tunnel splitting 
A appears because the level crossing is between states of the 
same fermion parity. 
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FIG. 4: Tunnel splitting at (^o = tt as a function of the in- 
duced charge. The dashed curve correspond to Eq. (|15[) . the 
solid curve to numerical calculations for the full Hamiltonian, 
for £j,2 ^5Ec = 25Em,i (with Em,2 = = Ej,i). 



The tunnel splitting A ensures that the energy of the 
ring evolves 27r-periodically if the flux $ is increased by a 
flux quantum h/2e in a time T$ which is long compared 
to Ta = ^i?M, i/A^. For r$ < Ta there is a significant 
probability exp(— T^/Ta) for a Landau-Zener transition 
through the gap, resulting in a 47r-periodic evolution of 
the energy. 

This limiting time scale Ta originating from quantum 
phase slips can be compared with the time scale Tp for 
quasiparticle poisoning. We require T$ small compared 
to both Ta and Tp to observe the 47r-periodic Joseph- 
son effect. For A > {hEM,i/TpY''^ one has Ta < Tp, 
so quantum phase slips govern. A recent experiment 
finds Tp ~ 2 ms in Al for temperatures below 160 mK.— 
Since Em,i will be well below 1 meV, one has Ta < Tp if 
quantum phase slips occur with a rate A//i higher than 
30 MHz. While quantum phase slip rates can vary over 
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many orders of magnitude due to the exponent in Eq. 
(fTS)) . typical values for a DC SQUID are in the GHz range. 

In conclusion, we have shown that Coulomb charg- 
ing effects do not spoil the 47r-periodic Josephson ef- 
fect in a superconducting ring, provided that all weak 
links contain Majorana fermions. Quantum phase slips 
at a weak link without Majorana fermions restore the 
27r-periodicity on time scales long compared to a time 
Ta, which may well be shorter than the time scale for 
quasiparticle poisoning. 

The origin of the protection of the Att periodicity if the 
entire ring is topologically nontrivial is conservation of 
fermion parityi^ (See Ref. |27j for a more general perspec- 
tive.) This protection breaks down if part of the ring 
is a trivial superconductor, because then the level cross- 
ing involves states of the same fermion parity and tunnel 



splitting by the charging energy is allowed (see Fig. [3]). 

We note in closing that the different stability of the 
47r-periodic Josephson effect in the two geometries of 
Fig. m examined here with respect to Coulomb charg- 
ing, extends to other parity-preserving perturbations of 
the Hamiltonian. For example, overlap of the wave func- 
tions of two Majorana bound states on the same island 
introduces a term -ffovcriap = *e7i72- For the lower panel 
of Fig. [I] this term leads to a tunnel splitting A = 2e 
which spoils the 47r-periodicity.— For the upper panel of 
Fig. [U A = because ^overlap preserves fermion parity. 
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